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Abstract. We revisit the question of the measured, unexpectedly large, heavy 
flavor suppression, RjiAiPT) li in nucleus-nucleus collisions at RHIC and 
compare two new theoretical approaches to the D- and _B-meson quenching. In 
the first model, radiative energy loss, coUisional energy loss and heavy quark- 
resonance interactions are combined to evaluate the drag and diffusion coefficients 
in the quark-gluon plasma and the mixed phase. These are applied in a relativistic 
Fokker-Planck equation to simulate the heavy c- and b-quark suppression rate 
and elliptic flow V2{pt)- In the second model, the fragmentation probability 
for heavy quarks and the medium-induced decay probability for heavy hadrons 
are derived. These are implemented in a set of coupled rate equations that 
describe the attenuation of the observable spectra from the coUisional dissociation 
of heavy mesons in the QGP. An improved description of the non-photonic electron 
Raa{pt) at RHIC can be obtained. In contrast to previous results, the latter 
approach predicts suppression of B-mesons comparable to that of D-mesons at 
transverse momenta as low as pT ~ 10 GeV. 



1. Introduction 

The detailed suppression pattern, RAAiPr), and elhptic flow, V2{pt)i of high- 
transverse-momentum hadrons is an important experimental signature of the quark- 
gluon plasma creation in heavy ion collisions [T]. Jet quenching for light mesons, such 
as TT, K and rj, at RHIC is well explained by radiative energy loss calculations [2j. It 
also gives the dominant contribution to the azimuthal asymmetry of hard probes [1] . 
In contrast, models [5] with a physically reasonable set of QGP temperatures and 
densities, predict a QCD heavy-quark energy loss which is too small compared to the 
measured suppression of single non-photonic electrons [31 13]. Therefore, it is critical 
to investigate new interaction mechanisms in the QGP that may be specific to heavy 

flavor [inmniin]. 



2. Heavy-flavor suppression in a combined transport + quenching 
approach 

Thermalization of heavy quarks in the QGP-heat bath has been recently studied in 
the framework of the parton-transport approach [6l El [8] . Large nuclear suppression 
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Figure 1. Left panel: preliminary results on the nuclear modification Raa{pt) 
for heavy c- and 6-quarks I10| from coUisional [5] and radiative energy loss [9] 
and quark-resonance interactions [g] [?]• For charm quarks, the PQCD AE 
contribution is shown separately. Right panel: elliptic flow V2{pt) for heavy 
c- and b-quarks for the same physics mechanisms I10| . 



and elliptic flow V2 result when employing a Fokker-Plank equation, 



df{p,t) _ d 



d_ 

dp J 



P,A{p,t) + —B,j{p,t) 



(1) 



dt dp' 

solved via an equivalent Langevin simulation. In Eq. ([T]) /(p, t) is the distribution of c- 
and 5-quarks and A{p, t) , B {p, t) are the drag / diffusion coefficients, respectively. It 
has been argued that strong coupling between the c- and 5-quarks and the QGP 
medium may be generated via quark-resonance interactions near the QCD-phase 
transition, T ~ Tc [5]. It is, therefore, important to study the interplay between such 
non-perturbative effects and the radiative [P] and coUisional [:5] heavy-quark energy 
loss. We evaluate [1(1 the contribution of these processes to the drag and diffusion 
coefficients. 



Aip,t) = ^^, B.,Mt) = -^-^^, (2) 
' pi St ' ' 2 6t ' ^ ' 



which are then applied in the relativistic Fokker-Planck equation. 

Results on the py-dependent suppression pattern of heavy quarks, -RaaIpt), are 
shown in the left panel of Fig. [1] Drag coefhcients are easily evaluated from the 
fractional momentum loss of heavy quarks, see Eq. The diffusion coefficients in 
this preliminary study were constrained from the fluctuation-dissipation relation. We 
observe that the suppression of charm quarks can be very large even in minimum-bias 
reactions of large nuclei and (charm) ^ (bottom). The high-px azimuthal 
asymmetry for minimum-bias Au-t-Au collisions is shown in the right panel of Fig. [T] 
We note that the generated V2 for 6-quarks is much smaller than that for c-quarks. 

One of the reasons for the large suppression in our current energy-loss 
implementation is that the Einstein fluctuation-dissipation relation induces minimal 
Gaussian fluctuations. These are significantly different from the ones in the 
probabilistic treatment of PQCD-energy loss [Tl[2l[5]. Future Langevin simulations of 
c- and 6-quark diffusion should include momentum fluctuations beyond the Einstein's 
relation and the decay of the heavy quark / hadron spectra into (e^ -I- e^) for direct 
comparison to the non-photonic electron observables at RHIC [TU]. 
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Figure 2. Left panel: suppressions of D- and _B-meson production via coUisional 
dissociation in the QGP in central Au+Au and Cu+Cu reactions at RHIC 
for gluon rapidity densities dN^/dy = 1175 and 350, respectively [2]. Right 
panel: quenching of the non-photonic electrons from the softened D- and B- 
meson spectra in central Au+Au collisions | 11| . Data is from PHENIX [3l I14| 
and STAR [1[12][13]. 

3. QGP-induced dissociation of heavy mesons 

In the perturbative QCD-factorization approach, the cause of the Umited single 
non-photonic electron quenching is identified as the small suppression of B-mesons, 
which dominate the high-pT e+e~ yields. Such models assume that the hard jet 
hadronizes in vacuum, having fully traversed the region of hot and dense nuclear 
matter, L^*^^ < 6 fm, and lost energy via radiative and coUisional processes [lj|2l|5]. In 
Ref. [11] we examined the validity of this assumption for different species of final-state 
partons and decay hadrons. For a, pr = W GeV pion at mid-rapidity Tform ~ 25 fm 
^ L^'^^ , consistent with the jet-quenching assumptions [Hll]. In contrast, B- and 
D- mesons of the same pt have formation times Tform ~ 0.4, 1.6 fm, respectively, 
<C L^*^^. Therefore, at the finite px range accessible at RHIC and LHC a conceptually 
different approach to the description of D- and B-meson quenching in A+A collisions 
is required, when compared to light hadrons. 

Motivated by this finding, in the framework of the GLV theory, we derive the 
coUisional dissociation probability of heavy mesons in the QGP [TT] : 

Pdixi^^O = [1 - Psixt^^O] > , Pdixf^^C = 0) = . (3) 
In Eq. ^ 2x/i^C — 2(/^^i/A)^ is the cumulative 2D transverse momentum squared 
per parton. The dissociation probability also depends on the detailed heavy meson 
light cone wave function. The dynamics of open heavy flavor production and 
modification in this model is represented by a set of coupled rate equations that 
describe the competition between b- and c-quark fragmentation and D- and _B-meson 
dissociation [TT] : 

dtf'^iPT, t) = -J^^^^IL^ + ^ /' dx \c^Q,H{x)f"{pT/x, t) , (4) 

dtf"iPT,t) = ^J^lifll±^+ / f\z^DH/Qiz)f'^iPT/z,t) . (5) 
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In Eqs. dm and © p{pT,t) = da'/dycPpr- For further details, see [II]. 

We solve this system of coupled rate equations numerically, using the same initial 
soft-gluon rapidity density dN^/dy as in the calculation of the 7r° quenching 2J in 
central Au+Au and Cu+Cu collisions at RHIC. Our results, including a study of a 
range of anticipated QGP densities at the LHC, are shown in the left panel of Fig. [2] 
Contrary to calculations that emphasize radiative and collisional heavy quark energy 
loss O [71 [To] , QGP-induced dissociation predicts i?- meson suppression comparable to 
or larger than that of £)- mesons at transverse momenta as low as pr ^ 10 GeV [llj . 
The heavy meson spectra are decayed into electrons (e++e^) using the PYTHIA event 
generator. Our results are shown in the right panel of Fig. [2] The predicted RaaIpt), 
which does not neglect the large i?-meson contribution, describes well the most recent 
heavy flavor quenching measurements at RHIC [12l [HI [T^ . We emphasize that such 
agreement between theory and experiment is not achieved at the cost of neglecting 
the contribution of the i?-mesons to the non-photonic e+e" spectra. 

4. Conclusions 

In these proceedings, we compared two new theoretical approaches [10( 111 1 to open 
heavy flavor modification in the QGP. Preliminary results on Langevin simulations 
of heavy quark diffusion, which include radiative energy loss, collisional energy loss 
and quark-resonance interactions, were shown. While an improved implementation of 
momentum fluctuations is required for quantitative comparison between data [Sj [4] 
and theory lOJ, we find normal suppression, i?^^(c) ^ and elliptic 

flow, U2(c) 3> V2{b), hierarchies as function of the heavy quark mass. Results 
on QGP-induced collisional dissociation of heavy mesons [H] were also shown. A 
good description dH [131 [H] of the large quenching of the inclusive non-photonic 
electrons |3l [4] is achieved by this model. A natural consequence of the approach 
developed in Ref. [11] is that B-mesons are attenuated as much as Z?-mesons at 
transverse momenta as low as pr ~ 10 GeV. We conclude that robust experimental 
determination of the dominant mechanism for in-medium modification of open heavy 
flavor would require direct and separate measurements of the B- and I?-meson Raa 
and V2 distributions versus pT and centrality in collisions of heavy nuclei. 
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